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ABSTRACT 

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large- 
scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared 
(FIR) [CII]158^m and [OI]63/im lines and mid-infrared H2 emission lines as tracers of gas cooling, 
and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, 
using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the 
ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating 
efficiency ([CII]158/im+[Oi]63^m)/PAH in the ring about 50% lower than in Enuc N and S. The 
average 11.3/7.7/im PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized 
PAHs, but no clear correlation with [CII]158/im/PAH(5.5 — 14/im) is found. PAHs in the ring are 
responsible for a factor of two more [Cll]158/im and [OI]63/im emission per unit mass than PAHs in 
the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can 
come from high intensity photodissociation regions (PDRs), in which case Go ~ 10 2 3 and nn ~ 10 3 ' 5 
cm~ 3 in the ring. For these values of Go and n# PDR models cannot reproduce the observed H2 
emission. Much of the the Hi emission in the starburst ring could come from warm regions in the 
diffuse ISM that are heated by turbulent dissipation or shocks. 

Subject headings: galaxies: starburst — galaxies: individual (NGC1097) — infrared: galaxies 
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1. INTRODUCTION 

The most important heating process in the neutral 
ISM is expected to be the photoelectric effect working 
on polycyclic aromatic hydrocarbons (PAHs) and small 
dust grains (iWatsonl Il972t iTielens fc Hoflenbachl 119851: 
IWolfire et all 119951: IHollenbach fe Tielensl I1999D . Inci- 
dent far-ultraviolet photons have energies high enough to 
eject electrons from dust grains {hv > 6 eV). They heat 
the gas via these photoelectrons, with a typical energy ef- 
ficiency of 0.1%— 1%. The efficiency is determined mostly 
by the ratio of UV radiation field to ga s density (Gp/n) 
which sets the charge of the gr ains (|Bakes fe Tielensl 
H994tlWeingartner fc Drainell200lD . The heated gas cools 
primarily via fine-structure lines such as [CII]158/im and 
[OI]63/im, with [CII]158/im being dominant at lower den- 
sities and temperatures. The photoelectric effect on dust 
thus couples the far-ultraviolet radiation field to the in- 
terstellar gas heating outside of H II regions. Such re- 
gions where the dust photoelectric effect dominate gas 
heating and/or chemistry are calle d photodissociation 
regions (jHollenbach fc Tielensl [1999) . Rotational tran- 
sition lines of Hi are also produced in photodissociation 
regions (PDRs), at the interface between the molecular 
gas and the atomic gas. The heating efficiency of the 
photoelectric effect is typically measured as the photo- 
electric yield per heating luminosity, using the flux ratio 
of the [CII]158/jm and/or [OI]63/im to the far-infrared 
(FIR) luminosity. A further discussion on the definition 
of the photoelectric heating efficiency can be found in 
Section l4~2l 

Observations of the FIR cooling lines in represen- 
tative samples of local star-forming galaxies were pi- 
oneered with the Kuiper Airborne Observatory, but 
were first studied in local s tar- fo r ming galaxie s with 
ISO (IMalhotraet all [19971 [20011 IHelou et al.l l200l 
iContursi et al.ll2002f) . These studies clearly showed that 
the flux ratio [CII]158/im/L(FIR) decreases by more than 
an order of magnitude (from 0.004 to less than 0.0004) 
for galaxies with high L(IR)/L{B) and/or warm dust 
temperatures while the [CII]158/im/PAH ratios show 
no such trends with dust temperature or luminosity. 
This means that the role of PAHs on gas heating de- 
creases with L(FIR)/ L(B), which could be attributed 
to PAH ionization, or that [CII] becomes a less im- 
portant cooling channel, resulting in more cooling via 
[OI]63/im. In "active" regions with high FIR luminos- 
ity, [OI]63/im becomes a more important co oling mecha- 
nism than [CII]158/mi (jMalhotra et al.l2001l) . The mech- 
anism responsible for heating these "active" regions is 
uncertain, but shock heating is one candidate, and en- 
ergy injection from deeply embedded star formation in 
dense regions is another. Models of the integrated emis- 
sion of [OI]63/im, [CII]158/im, Hi vibrational modes and 
other emission lines from Galactic and extragalactic star- 
forming regi ons have been used to diagnose these mech- 
amsms (e.g. IHollenbach fc McKed["l989l: Kaufman et~aTl 
2006) . Unfortunately the low spatial resolution of ISO at 
long wavelengths and the modest samples of sources de- 
tected in multiple lines make direct comparisons to mod- 
els extremely difficult. With Herschel-PACS and Spitzer- 
IRS, there is now the opportunity to use the PAH bands 
and the main infrared cooling lines to study the varia- 
tions of the photoelectric heating efficiency and diagnose 



the main heating mechanisms of the ISM in an environ- 
ment with a wide range o f ISM phases. 

With Herschel /PACS (|Poglitsch et al.l I2010D we are 
now able to target the most important cooling lines of the 
warm ISM on physical scales much smaller than previ- 
ously possible. The KINGFISH project (Key Insights on 
Nearby Galaxies: a Far- Infrared Survey with Herschel - 
PI: R. C. Kennicutt) is an open-time Herschel key pro- 
gram which aims to measure the heating and cooling of 
the gaseous and dust components of the ISM in a sample 
of 61 nearby galaxies with the PACS and SPIRE instru- 
ments. The far infrared (FIR) spectral range covered 
by PACS includes several of the most important cooling 
lines of the atomic and ionized gas, notably [CII] 158 fj,m, 
[OI] 63 fim, [OIII] 88 fj,m, [Nil] 122 fj,m, and [Nil] 205 
/Ltm. 

NGC 1097 is a Seyfert 1 galaxy with a bright star- 
burst ring with a diameter of 2 kpc and a strong large- 
scale stellar bar (tG^rineLalJ 119881: iKohno et all 120031: 
iHsieh et al.l I2008D with a length of 15 kpc. This bar 
may be responsible for driving gas into the central re- 
gion of the galaxy, forming a ring of gas near the Inner 
Lindblad Resona nce, triggering the fo rmation of massive 
star clusters (e.g. lAthanassou lal ll992T ). There is also gas 
inflow inside the ri ng possibly fueling the central super- 
massive black hole (jPrieto et all 120051: iFathi et al.l 120061: 
iDavies et al.l 120091) . triggering the formation of a com- 
pact star cluster seen near the nucleus. The star-forming 
ring and the nucleus of N GC 1097 are prominent in 
CO 1-0 (|Kohno et al.ll200l as well as ro- vibrational H2 
and hydrogen recombin ation lines (jReunanen et al.|[2002t 
iKotilainen et al.ll2000() . Tra cers of denser mole cular gas 
such as CO 2-1 and HCN (|Kohno et al.ll200l peak at 
the nucleus and where the dust lanes and the ring in- 
tersect. A relatively bright, central concentration of the 
HCN emission also coincides with the nucleus, caused ei- 
ther by relatively dense molecular gas (n ~ 10 5 cm~ 3 ) or 
by a strong X-ray radiati on field from the nucl eus that 
is affecting the chemistry (Mciicrink et al. 20071) . 

Far-infrared ph otometry of NGC 1097 using PACS was 
first presented in lSandstrom et all ()2010D . which showed 
that the ring dominates the FIR luminosity of the galaxy. 
Therefore, with PACS it is pos sible to isolate the ring 
contribution to FIR luminosity. iSandstrom et al.l ()2010h 
also showed that mid- and far-IR band ratios in the ring 
vary by less than ±20% azimuthally, indicating modest 
variation i n the radiation fi e ld hea ting the dust on 600 
pc scales. Sandstr om et al.l (|2010D provides a better es- 
timate of the total bolometric emission arising from the 
active galactic nucleus and its associated central star- 
burst. The properties of the d iffuse gas in t he cen tral 
ring have been investigated by IBeirao et al.l (|2010l ) us- 
ing the PACS spectrograph. They observed an enhance- 
ment of [OI] 63 ^m/[CII] 158 /im ratio and [Nil] 122 
/im/[NII] 205 fj,m in the ring, which could correspond to 
an enhancement in t he intensity o f radia tion field and 
density in the ring. iCroxall et al.l (|2012t ) also includes 
PACS observations of NGC 1097 to study general trends 
of the ISM heating and cooling. NGC 1097 is therefore 
an ideal system to study the physical conditions of the 
ISM in a nearby galaxy with both a starburst and an ac- 
tive nucleus. However, the FIR maps used in this study 
does not probe HII regions, nor the warm molecular gas, 
which can trace shock and AGN excitation, nor do they 
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TABLE 1 

Exposure times 



Region 


SL 


LL 


SH 


LH 




(5.5 - 14.5/zm) 


(14 - 38Atm) 


(10.5 - 19.5/W) 


(19 - 38/im) 


Nucleus & Ring 


2196 s 


2139 s 


6696 s 


930 s 


Enuc N 


3294 s 




5544 s 


2684 s 


Enuc S 


3294 s 




6804 s 


2684 s 



compare the Herschel data to the emission in the mid- 
infrared where the strong PAH features reside. These 
ISM phases are better studied using emission lines found 
in the mid- infrared (5 — 40/^m). 

In this paper we present a study of the properties of the 
heating and cooling processes of the ISM in the ring and 
bar of NGC 1097 based on spectral maps using PACS 
on board the Herschel Space Observatory and the In- 
frared Spectrograph (IRS) on board the Spitzer Space 
Telescope. These two instruments are necessary to cover 
the wavelength range from 5.5 — 230/im, where the most 
important tracers of gas heating and cooling are emitted. 
We will diagnose ISM properties such as PAH ionization 
and heating efficiency to answer the following questions: 
how does the gas heating efficiency vary with PAH ion- 
ization? Is there a variation of PAH heating efficiency 
between the ring and the regions at the bar? Is the me- 
chanical energy created by the inflow of gas into the ring 
an important component of the gas heating in the ring? 
In this paper we address these issues in the following 
sections: in §2 we present the observations and describe 
the reduction process of the Spitzer-IRS data into data 
cubes; in §3 we present the spectra and the spectral maps 
of the nucleus and both ends of the bar for several spec- 
tral features such as PAHs and H^ lines and use ratio 
maps to analyze the PAH emission distribution, the gas 
heating efficiency, and the cause of the warm Hi emis- 
sion; in §4 we discuss the implications of the results in §3. 
Our main conclusions are summarized in §5. Throughout 
the paper we use a distance of 19.1 Mpc for a cosmol- 
ogy Hp = 75 km/s/ Mpc, using the Tully-Fisher relation 
(|Willicket alJll997ft . 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 1097 was a Science Demonstration Phase tar- 
get observed with the Herschel — PACS instrument as 
part of the KINGFISH (Key Insights on Nearby Galax- 
ies: a Far Infrared Survey with Herschel; PI: Kenni- 
cutt) Open Time program. The galaxy was observed 
with the PACS integral field unit in both chop-nod and 
wavelength-switching modes. The observations and data 
reduction pr ocedure of the PACS spectroscopy data were 
described in IBeirao et a l. (2010). The observations use 
an 2 x 2 oversampled raster map strategy (step size 12"). 
The standard pipeline procedure for this, and part of the 
HIPE chop-nod pipeline, is to grid these four separate ob- 
servations onto a common Ra-Dec grid using a smaller 
pixel scale than the native 9". 4 x 9". 4 pixels size (in this 
case 3" x 3" sub-pixels). For updates on PACS reduc- 
tion and a discussion on the differences between obser- 
vations in chop-nod and wavelength-switching mode, see 
iCroxall et al.l (|2012[ ) , which also studies general trends of 
the ISM heating and cooling. 

NGC 1097 was observed with Spitzer-IRS in spec- 
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Fig. 1. — Footprints of the Spitzer/IRS spectral maps of NGC 
1097, overlaid on a IRAC Sfim image from the SINGS program 
(ID 159). The labels refer to the Ring+nucleus (blue), Enuc N 
(northern edge of the bar, in red) and Enuc S (southern edge of 
the bar, in green), analyzed in this paper. The Enuc N and Enuc 
S regions have only short-low observations, while the ring+nucleus 
have spatial coverage of the IRS with short-low and long-low IRS 
modules. 

tral mapping mode as part of SINGS (Spitzer Infrared 
Nearby Galaxy Survey; iKennicutt et al.l (|2003l )). under 
program ID 159. The central region which includes the 
ring and the nucleus was observed using all high- and 
low-resolution modules, and the edges of the large scale 
bar, which we refer to as the northern and southern ex- 
tranuclear regions (Enuc N and S), was observed using 
all high-resolution modules but with only SL (5 — 14/im) 
at low resolution. The central region was observed on 
July 18, 2004, and the Enuc N and S regions were ob- 
served on August 8, 2005. Table 1 summarizes all the 
observations. In Figure [1] we show all the footprints of 
the IRS observations performed on NGC 1097. We ex- 
tracted only the maps that contain both modules SL1 
and SL2 or LL1 and LL2. While the short-low (SL) 
module was oriented in a parallel direction relative to 
the bar, the long-low (LL) module was oriented in a 
perpendicular direction relative to the bar. The high- 
resolution modules, SH and LH were oriented at an angle 
of 45 deg relative to the bar. Spitzer-IRS spectra of the 
central region of NG C1097 has been published before by 
iBrandl et al.l ([2006D and the Spitzer-IRS observations of 
the nuclear regions f or all SINGS galaxies including NGC 
1097 are described in S mith et al.l (|2007bl ) (low resolution 
spectra) and lDale et al.l fj2009f ) (high-resolution spectra). 
However, these measurements were averaged over large 
areas. Here we publish the complete IRS spectral maps 
of NGC 1097 for the first time. The total exposure times 
are given in Table [T] 

We used CUBISM (|Smith et all I2007al) to build SL 
and LL spectral cubes from sets of mapping mode 
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BCD (Basic Calibrated Data) observations taken with 
Spitzer-IRS. In order to create continuum-subtracted 
maps of spectral lines an d PAH features, we use PAH- 
FIT (jSmith et al.1 l2007bl ). PAHFIT is a spectral fit- 
ting routine that decomposes IRS low-resolution spectra 
into broad PAH features, unresolved line emission, and 
continuum emission from dust grains. PAHFIT allows 
for the deblending of overlapping features, in particular 
the PAH emission, silicate absorption and fine-structure 
lines. This is not possible using CUBISM. Because it 
performs a simultaneous fit to the emission features and 
the underlying continuum, PAHFIT works best if it can 
fit combined SL and LL spectra across the full 5 — 40/im 
range for which it was intended. In order to provide 
a stable long wavelength continuum from which to esti- 
mate the true PAH emission, including the broad feature 
wings, we formed a combined SL + LL cube by rotating 
and re-gridding the LL cubes to match the orientation 
and pixel size of the SL cubes (1".85 x 1".85). By calcu- 
lating the flux in the new LL cubes using a bilinear inter- 
polation (in surface brightness units), we assure that the 
flux over the native LL pixel scale is conserved. Finally, 
a small average scale factor of 0.81, calculated by com- 
paring pixels in the spectral overlap region between SL 
and LL at approximately 14 microns, is applied to the SL 
data before running PAHFIT on the combined SL + LL 
cubes. The output of PAHFIT is saved for each feature, 
and is used to create the PAH feature maps at the native 
resolution of the SL data as shown in Fig. [3l All subse- 
quent comparisons of the PAH emission with the PACS 
data are derived from the SL and LL cubes smoothed to 
the resolution of the relevant PACS features, as described 
in the text, and not from the high-resolution PAH im- 
ages shown in Fig. [3] We also used CUBISM to build 
SH and LH cubes from sets of mapping mode BCD (Ba- 
sic Calibrated Data) observations and to make maps of 
[SIII]18.7/mi, [SIII]33.5/im, and [SiII]34.5/xm. 

We also used Spitzer-IRAC 8/xm and Spitzer-MlPS 
24/mr images of NGC1097, obs erved as a part of SINGS 
(for a complete description see lDale et al.1 (|2005j )).The 8 
fim map was correc ted for stellar light using a scaled 3.6 
fim image following iHelou et all (pOOl . 

3. RESULTS AND ANALYSIS 
3.1. Spectra 

In Figure[2]we present the average low-resolution spec- 
tra of the four regions highlighted in this study (ring 
(blue), Enuc N (red), and Enuc S (green)), and a small 
region centered in the nucleus (in gray). These spectra 
were produced by extracting a rectangular area from each 
spectral cube using CUBISM: 17" x 17" for the ring, and 
3.7" x 3.7" for the nucleus. For Enuc N, the areas are 
24.1" x 40.7" and for Enuc S, the areas are 16". 7 x 37". 
The spectrum of the ring is the only one that covers the 
full 5 — 38/im range sampled by the SL and LL data. The 
nuclear spectrum has been subtracted from the ring spec- 
trum in Figure [2]to avoid contamination from the central 
AGN. Our low resolution spectra are very similar and 
show the typical characteristics of star-forming regions: 
prominent PAH features, a possible absorption feature at 
9.7/im and a continuum that rises with wavelength. To 
make a quantitative comparison between these regions, 
we used PAHFIT to decompose the low resolution spec- 
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Fig. 2. — Spitzer/IRS low resolution spectra of the ring (blue), 
Enuc N (red), and Enuc S (green), and the nucleus (grey). All 
low resolution spectra include the short-low wavelength range (SL; 
5 — 14/im) but only the ring includes IRS long-low (LL; 14 — 38^tm). 
The Enuc N and Enuc S regions were not observed in LL. The 
labels indicate the positions of the main spectral features in the 
mid-infrared. 

tra. We measured the flux of the PAH features, ionic 
lines and H2 lines, and the results are listed in Tables 
??, and [3J One of the most prominent differences be- 
tween the nuclear region and the ring is the flux of the 
H2 lines at 9.1 fim and 12.3/xm compared to the PAH fea- 
tures. The flux ratio H 2 S(3)/PAH 7.7 fim is - 0.009 in 
the nucleus but only ~ 0.004 in the ring and ~ 0.001 in 
the Enuc regions. 

3.2. PAH emission distribution 

PAHs constitute between 10 — 20% of the inter- 
stellar carbon, and were first identified as a possi- 
ble carrier of the wide em ission features at 6.2 , 7.7, 
8.6, 11.3, and 12 7jim by iLeger fc Puget] (|1984D and 
lAllamandola et all (|1985l ). These bands are emitted by 
PAH molecules following photoexcitation. The same 
PAHs dominate the photoelectric heating of the dif- 
fuse ISM. Ionized PAHs are usually associated with 
the emission bands at 6.2, 7.7, 8.6/im, whereas neutral 
PAHs are responsible for the emission bands at 11.3, and 
12.7/ Ltm. The larger PA Hs e mit at longer wa velengths 
(e.g. lAllamandola et al.lll999HHonv et al.ll2001l ). For this 
reason the PAH ratios 11.3/7.7/xm can be used to diag- 
nose PAH ionization, and 6.2/7. 7/xm can be used to di- 
agnose grain size (e.g. iDraine fe Li 2001), although with 
caveats, as the 11.3/7.7/zm also varies slightly with size. 

In Figure [3] we present spectral maps of the continuum- 
subtracted total PAH emission between 5 — 14fim for the 
ring, the northern end of the bar (Enuc N) and the south- 
ern end of the bar (Enuc S). Each map was constructed 
by smoothing the 6.2, 7.7, 8.6, 11.3, and 12.7/im to the 
14/im resolution, assuming that the PSF at 14/im can be 
approximated to a gaussian with FWHM ~ 4".l. The 
map of the ring reveals that the average PAH surface 
brightness in the ring is a factor of 10 higher than in 
the nucleus and outer regions. The PAH surface bright- 
ness varies by a factor of ~ 4 throughout the ring, and 
peaks in a hotspot north of the nucleus. This is the same 
hotspot where a maximum of [OIl63/i m/[CII1158Mm was 
also observed by iBeirao et al.l ([2010), which was inter- 
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Fig. 3. — Maps of the total PAH emission between 5 — lAfim 
in the ring and nucleus (upper), Enuc N (center), and Enuc S 
(lower). Each image has a resolution matched to the IRS PSF at 
14/^m which we approximate as a 4".l FWHM Gaussian. The pixel 
scale of l".85/pixel. North is up. 

preted as a region with an intense radiation field. PAH 
flux in the nucleus is a factor of ~ f lower than in the 
ring, but as seen in Fig. 2 is still quite strong. We can 
compare this result to other studies which involve PAH 
observations in the nucleus on NGC f 097 at a higher spa- 
tial resolution. iMason et all (|2007t ) found no PAH 3.3um 
in the central 0".2, but they do find strong PAH f f .3um 
emission over a 3.6" region. As seen in Fig. 2 we also 
observe strong PAH emission in the nucleus over a simi- 
lar area, so our results are in agreement. In the Enuc N 
and S, the PAH surface brightness varies by a factor of 
~ 10, peaking near star forming regions. 

In Figure 2] we show spectral maps of the 11.3/7.7/xm 
PAH ratios of the ring and extranuclear regions at both 
ends of the bar. These maps were produced using 7.7/im 
and 11.3/im PAH maps smoothed to a 4".l resolution us- 
ing a gaussian kernel. The ratio increases up to a factor 
of 2 from a low of ~ 0.2 in the ring to a peak of ~ 0.4 
outside of the ring and in the nucleus. In the extranu- 
clear regions the ratio varies in the same proportion, with 
the lowest ratio coinciding with the star forming regions. 
However, in Enuc N the PAH surface brightness peak at 
the southern end of the map does not coincide spatially 
with a low value of 11.3/7.7/xm PAH ratio. This is sur- 
prising, but the reason is not clear, and indicates that 
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we approximate as a 4".l FWHM Gaussian. The pixel scale of 
l".85/pixel. North is up. 

the PAH ratios are not always a direct indicator of the 
starlight intensity. In Figure [5] we plot the 11.3/7.7/im 
PAH ratio as a function of 6.2/7.7/im using PAH 6.2, 
7.7, and 11.3 maps for the ring (blue), Enuc N (red), 
and Enuc S (green). These maps were convolved to the 
resolution at 11.3/im, 2". 7, assuming a gaussian PSF. 
Each point corresponds to a 3". 7 x 3". 7 region. The 
tracks indicate 6.2/7.7/im and 11.3 /7.7/ffli PAH ratio s 
for neutral and ionized PAH grains (|Draine 
and the dot-dashed line indicates where PAH grains have 
approximately 250 C atoms. The error bars indicate rep- 
resentative measurement uncertainties. The 11.3/7.7/im 
values for the ring are concentrated at ~ 0.23, but there 
is a wider dispersion for the Enucs. The 6.2/7.7/im ratio 
is on average 20% larger in the Enuc N and S than in the 
ring, but 20% lower in the nucleus than in the ring. This 
suggests that on average PAH grains may be smaller in 
the Enuc regions, but larger in the nucleus than in the 
ring. Table 2 shows the integrated fluxes of the main 
PAH features and PAH ratios in the ring, nucleus, Enuc 
N and Enuc S. 

3.3. Fine- Structure Lines and PAH Emission 

For a direct comparison between the PACS and the 
IRS data, we measured the flux of the [CII]158/zm and 
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Fig. 5. — Plot of the 11.3/im/7.7^tm PAH flux ratio as a func- 
tion of 6.2/im/7.7^tm PAH flux ratio. Each point corresponds to 
a 3". 7 X 3". 7 pixel region in the SL maps. The dashed lines indi- 
cate ionization states of the PAH molecules. Data from the ring 
is in blue, the nucleus is in violet, Enuc N is in red, and Enuc 
S is in green. The error bars indicate the average uncertainties 
for each region. The dash-dotted line indicates the position of the 
molecules formed b y 250 carbon atoms. This figure is adapted from 
IDraine h O <|200lT ). 
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Fig. 6. — The 12" X 12" square apertures selected for measuring 
the [CII]158/im, [OI]63/jm, PAH (5.5 - 14/W), 70/mi, and 100/im 
fluxes in the ring (blue), Enuc N (red), and Enuc S (green), overlaid 
on a IRAC 8/rni image of NGC 1097. North is up. 

[OI]63/im fine-structure lines, the total flux of the PAH 
emission features between 5.5 — 14/xm, the 11.3/7.7/xm 
PAH ratio, and the 70/im/100/xm ratio in the ring and 
extranuclear positions over a grid composed of 12" x 12" 
square regions, shown in Figure [HI We chose these regions 
to properly sample the flux at the longest wavelength 
where the FWHM is close to 11". Before measuring the 
total PAH flux and the 11.3/7.7/im PAH ratio, we made 
convolution kernels to smooth each PAH feature map 
at 6.2,7.7,8.6,11.3,12.6, and 14.1/im to the Herschel 
beam at 160/xm, using PSFs calculated bv lAniano et al.l 
(|2011f) . and theoretical PSFs for IRS made with STiny- 
Tim for each wavelength above. We then rebinned all 
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Fig. 7.— Top: Plot of the ratio [CII]158Mm/PAH (5.5 - 14/^m) 
as a function of the 70/xm/100/im flux ratio. Bottom: Plot of the 
ratio [CII]158/im/PAH (5.5 - 14/im) as a function of the 11.3/7.7 
PAH ratio. Each point was calculated for a resolution of 11" per 
beam. 
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Fig. 8.— Plot of the ratio ([CII]158/im+[OI]63^m)/PAH (5.5 - 
14/^m) as a function of the 70/mi/100/^m flux ratio. Each point 
was calculated for a resolution of 11" per beam. 

maps to match the pixel size used in the PACS FIR spec- 
tral maps (3"/pixel). The PAH 7.7/zm and the 11.3/xm 
fluxes were measured summing all pixels within each re- 
gion in Fig. |51 We also used [OI]63/^m maps and the 
70/zm and 100/Ltm PACS images convolved to the PACS 
160/zm PSF using the convolution recipe described in 
lAniano et ail 1)20111) . Due to low S/N, [OI]63Aim could 
only be measured in the ring and Enuc S. The total 
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luminosity of [CII]158/im and [OI]63/im in the ring is 
L[cii] ~ 1-4 x 1O 8 L and i [0/ ] ~ 5.4 x 1O 7 L , re- 
spectively. In Table 0] we present all the [CII] , [01] , and 
PAH fluxes measured for all the regions, as well as ionic 
lines measured at hi-res, [SIII]18/im, [SIII]33.5/tm, and 
[SiII]34.5/mi. 

In Figure [7] (upper) we plot the ratio [CII]158/im/PAH 
(5.5 -14/im) as a function of 70/im/lOO/tm flux ratios, 
based on PACS photometry for the ring and the two ex- 
tranuclear maps at a resolution of 11" per beam. The 
70/im/100/im flux ratio is a commonly used indicator 
of dust temperature. Each datapoint is one of the 12" 
square regions (~ 1.1 kpc). There is a clear differ- 
ence between the ring and the extranuclear positions in 
both [CII]158/mi/PAH (5.5 -14/im) and 70/im/100/mi. 
We performed a KS test to determine the significance 
of this difference. The D statistic is 0.85, with a 
probability of 0.029% of the ring and Enuc S belong- 
ing to the same population. The difference between 
the two regions is therefore significant. For the ring, 
the median [CII]158/xm/PAH (5.5 -14/im) is ~ 0.034 
and for the Enucs ~ 0.053, about a factor of ~ 1.6 
above the median ring value. The 70/iin/100/im ratio 
ranges between 0.35-0.55 for Enuc N and S and be- 
tween 0.70-0.80 for the ring. The dust is warmer and 
the [CII]158/im/PAH (5.5 -14/im) is lower in the ring 
than in the Enucs. If we consider the total infrared 
emission between 5 — 10 /im, F (5.5 -10/im ), we can 
compare our results with iHelou et al.1 (|200lD . in which 
[CII]158/tm/F (5.5 -lQ/im) is constant, suggesting that 
photoelectrons liberated from PAH grains were responsi- 
ble for heating the gas. Our values for [CII]/F(5 — 10/im) 
are ~ 50% lower in the ring than in the Enucs. In Fig- 
ure [7] (lower) we plot the same ratio [CII]158/tm/PAH 
(5.5 -14/im) now as a function of the 11.3/7.7/im PAH 
ratio, at a common resolution of 11" per beam. The 
distribution of 11.3/7.7/im PAH ratios in the ring is nar- 
rower than in the Enuc regions, and has a median value 
of 0.25, lower than the median value of 11.3/7.7/im in 
the Enuc regions, which is 0.28. This could be due to 
the fact that the PAH emission in the ring is more ion- 
ized. We performed a Spearman rank correlation test, 
and the result is 0.33, revealing a weak correlation be- 
tween [CII]158/im/PAH (5.5 -14/xm) flux ratio and the 
11.3/7.7/im PAH ratio. However, this apparent corre- 
lation is only due to the lower [CII]158/im/PAH (5.5 - 
14/im) values in the ring, whereas the Enuc values of 
[CII]158/mi/PAH (5.5-14/mi) appear to be insensitive to 
the local 11.3/7.7/im flux ratio. Therefore it is not clear 
that the difference in [CII]158/tm/PAH (5.5 -14/im) be- 
tween the ring and the Enucs is due to a variation of 
the PAH ionization. Note that the 11.3/7.7/im data- 
points shown in Figure [7J were measured over a much 
wider area, and therefore are an average of the data- 
points in Figures 4 and 5. This accounts for the narrower 
range of 11.3/7.7/tm values in Figure [71 [OI]63/tm can 
also be an important coolant in warmer regions. In Fig- 
ure [8] we plot the sum of the [CII]158/mi and [OI]63/im 
emission over the PAH flux between 5 — 14/tm against 
the 70/im/100/tm flux ratio. We measured the [OI]63/tm 
fluxes after smoothing the [OI]63/im image to the 160/im 
PACS PSF (a resolution of app roximately 11"), using 
kernels from lAniano et ahl (|2011[ ). Here we exclude Enuc 



N, due to low S/N in the [OI]63/tm line. The gap between 
the Enuc S and the ring in ([CII]158/im+[OI]63/im) /PAH 
is smaller than in [CII]158/mi/PAH, with a median of 
~ 0.047 for the ring and - 0.070 for Enuc S. [OI]63/im 
has therefore a much larger role in the cooling of the gas 
in the ring than in Enuc S. This is consistent with the 
behavior of ( [CII] 1 58/im+ [Oil 63 / im) /P AH in the galax- 
ies as reported by IHelou et al.l (|2001[ ) . Both the ring 
and the Enucs in NG C 1097 are within the range of 
parameters studied in Hcl ou et al.1 (|2001l ). They cal- 
culate [CII]/F(5 — 10/im) for a sample of star- forming 
galaxies using global averages, and shows no trend with 
60/im/100/tm. However, they find a large scatter in 
( [CII] + [OI] )/F(5 - 10/tm) as a function of 60/mi/100/iin. 
What we are uncovering is the inherent variations within 
galaxies that might drive the scatter in this relation. We 
should note that all the measurements are affected by 
the uncertainty on the PACS beam si z e and the shape of 
the PSFs, quantified in lAniano et al.l (|2011l) . 

Another important cooling line is [Sill] at 34.8/zm. 
With a critical density ~ 10 5 cm" 3 for collisions with 
neutral hydrogen and ~ 10 2 cm -3 for collisions with elec- 
trons, and a low ionization potential of 8.15 eV, this line 
comes from a wide variety of regions, such as ionized gas 
and warm atomic gas, photodissociation regions a nd X- 
ray dominated regions (jHollcnba ch fe Tie lcns 1999). HII 
regions represent a large fractio n of [SiII]34.8 /xm e missio n 
in the ring (between 30%-60%. lNagao et al.l (e.g. 12011ft ). 
We used an [SiII]34.8/tm image of the ring and Enuc 
S, and convolved it to the 160/im Herschel-PACS beam 
using theoretical Spitzer — IRS PSFs mad e with STiny- 
Tim and Herschel PSFs calculated by lAniano et al.l 
(pOll . From these images we measured a luminosity 
L[sni] ~ 4 x 10 8 i Q in the ring and L\ SiII } ~4x 10 7 L Q 
in Enuc S, so L^Sill] ~ 1.6 x iron- This extra cooling 
would raise the ring and Enuc S in Figure |SJ 

In Figure [5] we show a plot of [CII]158/im/FIR as a 
function of the 70/im/100/tm flux ratio, with the dat- 
apoints calculated at a common beam size (11"). FIR 
is the far infrared luminosity, estimated as FIR=TIR/2, 
with TIR as the total infrared luminosity, wh i ch wa s cal- 
culated using the formula by IDale fc Heloul (|2002t ). re- 
placing the 70/iin, and 160/tm flux by PACS 70/im and 
160/im fluxes. The ring has a [CII]158/im/FIR which 
is a factor of ~ 3 lower than the Enuc regions, and 
this difference is much larger than the factor of ~ 1.5 
in [CII]158/tm/PAH (5.5 -14/im) between the ring and 
the Enucs. This effect has been associated with re- 
gions with high Go/n wher e other cooling lines are en- 
hanced, such as [OI]63/tm ( Malhotra et al.ll200lD . The 
[CII]158/im/FIR in the ring is similar to the value in 
other star-forming galaxies observed by ISO such as NGC 
1482, which has [CII]158/xm/FIR ~ 0.004 and NGC1569 , 
with [CII]158/tm/FIR~ 0.0027 (|Malhotra et al.ll2^0l . 
In Figure [9] (bottom) we plot [OI]63/tm/FIR for the ring 
and Enuc S. We see that [OI]63/im/FIR is still lower for 
the ring, but [OI]63/iin contributes at most ~ 25% of 
the [CII]158/im+[OI]63/im flux in the Enucs and - 30% 
of the [CII]158/im+[OI]63/tm flux in the ring. Although 
some of the cooling indeed goes through [OI]63/im, it 
cannot fully compensate for the difference between the 
ring and the Enuc S in [CII]158/im/FIR. The decrease 
in [CII]158/im/FIR can be caused not only by cooling 
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Fig. 9.— Top: Plot of [CII]158/im/FIR as a function of the 
70/^m/100/im flux ratio. Bottom: Plot of [OI]63/im/FIR as a func- 
tion of the 70^tm/100/im flux ratio. 

via other emission lines, but also by absorption of ion- 
izing photons by dust, low PAH abundance due to PAH 
destruction or coagulation into larger grains, collisional 
suppression of [CII]158/xm, or the [CII]158^m being op- 
tically thick. 

3.4. Warm Hi emission 

Figure [TU] shows flux maps of the Hi rotational lines 
from S(0) to S(3) at their native resolution, extracted 
pixel-by-pixel with PAHFIT. Most of the S(0) flux is 
emitted from two spots in the ring that coincide with 
the intersection of the bar with the ring. The S(l) line 
flux is also emitted from these spots, but also peaks at 
the nucleus. In S(2) and S(3) the flux peaks at a region 
in the NE part of the ring, where the main Ha and FIR 
lines are emitted ()Beirao et alj l2010h . but there is also 
some emission from the nucleus. 

With the H2 emission maps, we can estimate the 
gas temperature, column density, and mass distributions 
by fitting the Hi fluxes using the Boltzmann equation 
( Rigo poulou et al] 120021) . assuming that the ortho/para 
ratio is in local thermodynamical equilibrium (LTE). In 
Figure[TT]we plot the excitation diagrams for the ring and 
nucleus of NGC 1097 (upper), and three selected (~ 10") 
regions on the ring (lower), with fits to the warm and hot 
Hi phases. In this fit, 100 K was set as the lower temper- 
ature limit, as this is approximately the lowest temper- 
ature in a collisionally excited Hi gas. Before the fit we 
convolved the S(3), S(2) and S(l) images to a gaussian 



with FWHM = 7", which is the resolution at 28/im. In 
Table [5] we present the properties of the warm Hi gas 
in the ring, nucleus, and three 10" diameter regions in 
the ring indicated in Figure 1101 The lower temperature 
component in the nucleus has ~ 2% of the total Hi lu- 
minosity of the nucleus, and therefore almost all the Hi 
lines can be fitted with a single temperature. The ring 
is fitted with two components, a cool component with 
a temperature of 119±11 K and a hot component with 
467±79 K. In the ring, the cool component emits 20% of 
the Hi flux, so there is a large quantity of gas in the ring 
at ~ 100 K that is not present in the nucleus. The lower 
plot in Figure [TT] shows that the hot component varies 
between T ~ 375 and 470 K. Of the three selected re- 
gions, Region 2 has the highest fraction of emission from 
the cool component, ~ 16%, and the lowest temperature 
of the warm component, perhaps indicative of a reservoir 
of colder gas in this region. 

4. DISCUSSION 

4.1. Physical conditions in the gas 

To explore the physical condi tions of the gas in the 
ring and Enuc regions we use Kaufman et all {2006) 
PDR models to determine the intensity of the radia- 
tion field and the density of the neutral gas using the 
[OI]63/im and [CII]158/im fine structure lines. Our mea- 
surements of [CII]158/im include an ionized gas contri- 
bution, where the gas is emitted from diffuse ionized 
regions. Therefore one needs to subtract the ionized 
contribution to [CII]158/im. This contribution can be 
estimated from measurements of the [Nil] 122/im and 
[NII]205/xm lines, which arise in diffuse ionized gas and 
HII regio ns. The [NI I ]122/u m/ [Nil] 205/im values pub- 
lished by iBeirao et all |2010l) were affected by an incor- 
rect Relative Spectral Response Function (RSRF) esti- 
mation, which affects the calibration of the continuum 
at A > 200/im. To correct this effect we have to multiply 
the [NII]205/im flux by a correction factor of 4.5, which 
results in a [NII]205/xm/122/Ltm~ 0.25 in the ring. 

Using the estimated [NII]205/xm/[NTI]122/mi value for 
the ring, we first calculate the fraction of [NII]122/Ltm 
from diffuse ionized gas (n e = 10 cm -3 ). iRubin et al] 
(fl99l calculate [NII]205/im/[NII]122/xm~ 1 for n e = 10 
cm" 3 and [NII]205/xm/[NII]122/im - 0.1 for n e = 10 3 
cm -3 . The average value of the [NII]205/im/122//m in 
the ring is ~ 0.25. Therefore we estimate that the diffuse 
medium contributes ~ 17% of the observed [Nil] 122/im 
emission in t h e ring . Using C + a n d N + abundances from 
| Mever et al] (119971). |Sofia et all (|2004D . A values from 



Galavis et al.l (fl997l. Il998f) and collision streng ths from 
iHudsonfc B ell (2005); Blu m k Pradhanl (|1992fi . we de- 
rived a [CII]158/im/[NII]122/mi~ 4.8 in diffuse ionized 
gas and [C/I] 158/ [JV7I] 122 - 0.8 in the dense ionized 
gas. If 17% of the observed [NII]122/xm is emitted from 
the diffuse ionized gas, [CII]158/xnij „ ~ 1.5[NII]122/im. 
We can use this estimate to calc ulate the ioni z ed ga s 
contribution to [Gil] 158jLim. From [Beirao et all (2010). 
[CII]158/xm/[NII]122~ 4.5 in the ring, which suggests 
that ~ 33% of the [CII]158/im flux in the ring is emit- 
ted from ionized gas and ~ 67% from PDRs. Note 
that this estimate is affected by the uncertainty on the 
[NII]205//m/122/iin value, which is ~ 50%. 
Similarly, we can determine the maximum ionized gas 
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Fig. 10.— Map of the H 2 S (0) and H 2 S(l), H 2 S(2), and H 2 S(3) at 28.2/^m, I7.0fim, 12.3£tm, and 9.7^m respectively. "Region 1", 
"Region 2", and "Region 3" denote three regions selected along the ring (see description in text). All images are at their native resolutions, 
i.e. 6". 8, 4".l, 3".l, and 2". 4, respectively if we approximate the IRS PSFs as Gaussians. North is up. 



fraction of the [CII]158/im emission in Enuc S, and com- 
pare to the value for the ring estimated above. Because of 
low S/N, no measurement of [NII]205^m was possible in 
Enuc N and S and no [NII]122/im was possible in Enuc N. 
However, we can use [SIII]18/im/33/im ratio as a tracer 
of density for regions with n e > 100 cm~ 3 , and use the 
density estimated from this ratio to esti mate the cor- 
respo nding [NII]122^m/[NII]205/xm (see iCroxall et all 
Mm . We measure [SIII] 18/33- 1.2 in Enuc S. This 
corresponds to [NII]122/im/[NII]205/im— 8. Assuming 
10 cm -3 for the diffuse component and 10 3 cm -3 for the 
dense component, we can estimate the diffuse gas contri- 
bution to [NII]122/rni/[NII]205^m in Enuc S as ~ 2.8%. 
Using these values, we get [CII]158/zm/[NII]122~ 0.9 for 
the ionized component in Enuc S, which is ~ 5% of the 
obser ved ratio, [CII]158/im/[NII]122~ 17 (jBeirao et al.1 
2010). Therefore, the ionized gas component accounts 
for ~ 5% of the [CII]158/im flux in Enuc S. Enuc S there- 
fore seems to have a much lower fraction of the observed 
[CII]158/im emission coming from ionized regions than 
does the ring. 

Using the neutral gas component of the [CII]158/im 
flux and the [OI]63/im flux for the ring and Enuc S, we 
estimate Go and i%h throughout these regions. We use 
the notation th at Gn is meas ured in units of the Habing 
radiation field (Habing: JL968J) . We first assume that, af- 
ter subtracting the ionized component, all the remain- 
ing [CII]158/^m emission and all [OI]63/^m and FIR flux 
is emitted from PDRs irradiated by m assive stellar clus- 
ters. Following IKaufman et al.l (|2006l ). the grid in Figure 



Q2] shows a variation of the incident radiation field Go be- 
tween ~ 80 in the Enuc S region and 800 in the ring. The 
gas density is 10 3-3 - 5 cm~ 3 in both regions. These val- 
ues are similar to local star- forming spiral galaxie s such 
as NGC 695 and NGC 3620 (jMalhotra et al.|[200l . The 
density is a factor of 10 lower than in M8 2, but the value 
of Gn in the ring is comparable to M82 (jKaufman et al.l 
1999). The larger Go in the ring presumably signals 
an increased ionization parameter, which measures the 
number of ionizing photons per hydrogen atom. This is 
consistent with our estimate above that — 33% of the 
[CII]158^m originates from ionized regions. 

However, when considering detailed fits of the dust 
SED of the ring instead of just the FIR emission line 
ratios, a different picture of the intens i ty of the ion- 
izing field emerges. Using iDraine fc Lil (|2007f ) models, 
lAniano et ail (|2011D estimate that only 25% of the TIR 
is being emitted from regions where Go > 10 2 . For Enuc 
S this fraction is only about 8%. Therefore, if we use 
these estimates for the infrared em i ssion together with 
the PDR models of [Kaufman et ail (|2006[ ). we predict a 
lower ionizing radiation field intensity for the ring, much 
closer to a value of Go — 10 2 3 , with a corresponding 
uh — lO 3 5 cm -3 . These significantly different estimates 
for the Go and nu in the ring PDRs have important im- 
plications for the production of the warm H2 emission 
we measure in the IRS data, and we will carry both of 
these estimates forward when discussing the source of 
the H 2 emission in Section 4.3. For Enuc S lAniano et al.l 
(120111) predict a PDR fraction of the TIR flux of ~ 8% in 
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Fig. 11. — Top: mid-infrared H2 excitation diagrams for the 
nucleus and the ring of NGC 1097. Bottom: excitation diagrams 
for the th ree regions selected along the ring and highlighted in 
Figure 1101 The lines indicate the best fit to the H2 datapoints. 
(see Table 0. 

Enuc S. This will make (F(\Cll])+F(\O l]))/FIR > 0.1 
in Enu c S, out of the range modeled by Kaufman et all 
(2006) for PDR regions. All the line emission could be 
coming from the diffuse ISM and not PDRs in Enuc S 
and therefore we cannot predict the Go and nu- 

Recently using Herschel/P ACS to observe a number o f 
infrared starburst galaxies, IGracia-Carpio et all (|201lD 
established a correlation between [CII]158/im/FIR and 
regions with high star formation efficiency, Lfir/ Mjj 2 , 
where [CII]158/xm/FIR drops only when Lfir/Mh 2 > 
SQLn jMp-,. Based on PACS 70/xm and lOO/itm fluxes 
from ISandstr om et al.l (|20100 and CO 2-1 fluxes from 
IHsieh et all (pOOl (assuming X CO = 3.0 x 10 20 cm 2 (K 
km s 1 ) 1 , and R21 ~ 1.3), we calculate Lfir/Mh 2 ~ 
50Lq/M q in the ring of NGC 1097, lower than the value 
of Lfir/Mh 2 for which [C II]158Atm/FIR decre a ses fo r 
the luminous starbursts in IGracia-Carpio et al.l (|2011l ). 
Although [CII]158/im/FIR is lower in the ring than 
Enucs, it is consistent with other starburst galaxies hav- 
ing similar values of Lfir/Mh 2 - 

4.2. Gas heating efficiency and PAH abundance 

Theoretical studies have shown that the extent of 
grain charging, and therefore the efficiency of photo- 
electric heating, is not only dependent on environmen- 



FlG. 12.— Plot of [CII]158Atm/[OI]63Atm as a function of 
([CII]158/xm+[OI]63/xm)/FIR luminosity ratios, emitted from a 
PDR according to a gri d of Go and njj, for the ring and Enuc 
S, as defined in Section 13.31 For clouds illuminated on all sides 
by Go, the observed ([O I]63^m + [C II]158/im)/FIR ratio is the 
ratio of the escaping intensities of line photons to grain continuum 
for each cloud. This corresponds to I([0 I] 63 /im + [C II] 158 
fim)/2 X Ifir, where Ipjn = 2 X 1.3 X 10 _4 Go ergs cm -2 s _1 
sr _1 . Only [OI]63/^m data from Enuc S and the ring are available. 
The error bars represent the average uncertainties derived from 
photometric errors. 

tal conditions, b u t also on grain size and grain species. 
iBakes fc Tielensl (|1994t ) estimated that approximately 
half of the heating is from grains with less than about 
1500 C-atoms (~ 15l). Using IRAS and ISO ob- 
servations of the [ CII1158^m, [OI]63/im and H 2 lines, 
lHabart et al.l ([2004) quantified the photoelectric efficien- 
cies attributed to the PAH, VSG (very small grains) 
and BG (big grains) populations in the p Oph cloud, 
finding cpah = 3%, eysc = 1%, an d £bg = 0.1% 
respectively. This suggests that PAHs are the most 
efficient gas heating channel in star-forming galaxies 
and t herefore drive t he heating of the neutral g as re- 
gions (IWatsonl 119721: IHollenbach fc Tielensl |1999f ). In- 
deed, IHelou et al.l (|200lh showed that the PAH contri- 
bution to the gas heating appears constant in galax- 
ies with increasing i 7 '(60/im)/i 7 '(100/im), and dominant 
over the contribution of large grains, based on the rela- 
tively constant [CII]158/im/PAH against a rapidly falling 
[CII]158^m/FIR. This suggests that PAHs dominate the 
heating of the [CII] -emitting gas (via the photoelectric ef- 
fect), wh ile large grains domin ate the FIR emission from 
galaxies ()Malhotra et al.ll2001ft . 

A decrease in the [CII]158/im relative to the FIR emis- 
sion could also be caused by a decrease in the photoelec- 
tric gas heating efficiency associated with highly charged 
grains. While the [CII]158^m/PAH (5.5-14/mi) is lower 
in the ring than in the Enucs in NCC 1097, we find no 
correlation between [CII]158^m/PAH (5.5 — 14/im) and 
the 11.3/7.7/im PAH ratio, which measures PAH ion- 
ization. Although the average 11.3/7.7/Ltm PAH ratio 
is lower on average in the ring than in the Enucs, the 
regions in the Enucs with similar 11.3/7.7/mi PAH ra- 
tios have clearly higher [CII]158/im/PAH (5.5 — 14/im). 
This means that the variation of [CII]158/im/PAH 
(5.5 — 14/im) between the Enucs and the ring can- 
not be attributed solely to PAH ionization. How- 
ever, we have seen in Figure [8] that [OI]63/im par- 
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tially compensates for the decrease of [CII]158/im/PAH 
(5.5 — 14/xm) ratio in the ring. As seen in Section 3.3, 
([CH]158/xm+[Oi]63/xm)/PAH in Enuc S is - 1.5 times 
higher than in the ring. After subtracting the ionized 
gas contribution to the [CII]158/xm emission, the heat- 
ing efficiency ([CH]158/im+[OI]63/im)/PAH in Enuc S is 
~ 1.9 times higher than in the ring. So by subtract- 
ing the ionized gas component, we obtain an estimate of 
the difference in the gas heating efficiency between the 
ring and Enuc S which is much larger than our previous 
estimates. 

Since the [SiII]34.8/im line is also a major coolant 
we can include this emission in the overall ener g y bud - 
get of NGC 1097. Following Kaufman et all il2006h . 
the [SiII]34.8/xm flux emitted from PDRs increases rel- 
ative to the [SiII]34.8/im flux emitted from HII re- 
gions as the electron density n e increases. The 
[SiII]34.8/im line can be as luminous as the [OI]63/xm, 
and ([CII] + [OI] + [SiII])/PAH(5.5 -14/xm) is similar in 
the ring and in Enuc S. However, this ratio would 
not represent accurately the gas heating efficiency, as 
it includes also the contributions to [CII]158/im and 
[SiII]34.8/xm from the ionized gas. Therefore we can 
only estimate upper limits on the photoelectric heating 
efficiency of the PAHs in the PDRs using ([CII]„ eut + 
[01] + [SiII])/PAH(5.5 - 14/xm) for the ring and Enuc 



S, respectively, with [CII 
gas component of the [CII 



neu t representing the neutral 
158/im emission, calculated in 
Section 14.11 This results in a PAH photoelectric heat- 
ing efficiency of 5.2 ± 0.3% in the ring and 7.7 ± 0.8%. 
Thus we have seen that adding [SiII]34.8/xm brings down 
the difference in photoelectric heating efficiency between 
the ring and Enuc S to 3a for ([CII]„ eut + [01] + 
[Sill] )/PAH(5. 5 - 14/xm). 

Using IRAC 8/xm and MIPS 24 /xm photometry o f 
the Large Magellanic Cloud (EMC), [Rubin et al.l ([20090 
noted that the ratio between the [CII]158/xm and 
the total infrared flux TIR decreases with decreas- 
ing F(8/xm)/F(24/im), which they use as a measure of 
PAH/VSG flux in regions with bright PAH emission. 
Variations of the F l/ (8/j,m)/F !y (24^m) ratio arise mainly 
from the destruction of PAHs in HII regions and heat- 
ing of VSGs to high temperatures in regions with in- 
tense radiation fields. If PAHs are destroyed, their emis- 
sion in the infrared (at 8/xm) decreases relative to the 
infrared emission from other grains (at 24/im). In re- 
gions with Go > 100, F v (24/xm) /Ftir will increase with 
starlight intensity, and the emission at 24/im will in- 
crease relative to 8/im. To verify if this is the case for 
NGC 1097, in Figure [TBI (upper - ) we present a plot of the 
([CII] + [OI])/PAH(5.5 - 14/xm) ratio for the ring and the 
Enuc N and S as a function of the F„(8/jim.)/ (24/im) 
ratio, calculated from the IRAC 8/im and MIPS 
24/xm maps convolv ed to the PACS 160/xm beam size 
(jAniano et alj 120111 ). The average F^S/xmVFV (24/xm) 
ratio in the ring in NGC 1097 is - 0.37 ± 0.015, the 
same found for locations in the LMC with the high- 
est [CII] 158/im surface brightness, whereas the average 
values of F„(8/im)/F„(24/im) for Enuc N and S are ~ 
1.1 ±0.1, which correspond to the average value observed 
in the LMC for similar [CII] 158/im surface brightness. 
The ring has both lower ([CII]158/xm+[OI]63/xm)/PAH 
and F J/ (8/im)/F^(24/(m) than the extranuclear regions. 
Therefore, the low ([CII] neut +[OI])/PAH(5.5 - 14/xm) 



and F v {d>iMn)/F v (24/xm) might be due both to a lack 
of PAH grains, and higher starlight intensity. We 
should note that F u (8fim)/ F^ (24/im) is not an inde- 
pendent parameter, as it is strongly correlated with 
F„(60it m)/i^(100/im). 

Using [Draine fe O ([2007Q models we can calculate the 
PAH mass fraction, defined as the fraction of grains 
with 1000 C atoms or less to all dust grains, as a re- 
sult of calculating infrared emission spectra for dust 
heated by starlight. In Figure [13] (lower) we present 
a plot of the ([CII] ne „ t + [OI])/PAH(5.5 - 14/xm) ratio 
for the ring and the Enuc N and S as a function of 
the PAH mass fraction (JpaHj as measured in the (/pah 
maps of NGC1097 constructed by lAniano et al.l ([2012). 
The trend for ([CII]„ e „ t +[OI])/PAH(5.5 - 14/xm) as a 
function of (/pah is similar to the trend seen for the 
F lJ (8fim)/F l , (24/im) ratio, with (/pah ~ 0.018 in the 
ring and ~ 0.035 in the Enuc regions, a difference of 
a factor of 2. The median (/pah in the SINGS sample 
is also 0.035 (jDraine et al.ll2007t) . This suggests that 
the variation of the F v {%iJ,m)/F v (24/xm) between the 
ring and the Enucs in NGC 1097 is a consequence of 
the variation of the PAH abundance relative to other 
types of grains. As long as PAHs dominate photoelectric 
heating and therefore [CII]158/xm and [OI]63/tm cool- 
ing, the ([CII]„ eut + [OI])/PAH(5.5 - 14/xm) ratio should 
not change with PAH abundance. As seen in Figure 
[H QCII] ne „t+[OI])/PAH(5.5-14/xm) decreases with de- 
creasing (XpAH and F u (S/inij/F^ (24/xm) ratio. This could 
be a consequence of an increase of the role of VSGs on 
the photoelectric gas heating, as the PAH mass fraction 
decreases and the VSG abundance increases faster than 
the PAH abundance. A significant portion of [CII] 158/im 
and [OI]63/im could be heated by VSGs, increasing 
([CII]„ eut +[OI])/PAH(5.5 - 14/xm). On the other hand, 
it may be that when (/pah increases, the PAH size dis- 
tribution shifts to smaller sizes, which may have higher 
photoelectric heating efficiencies. 

Figure [131 can also be interpreted by considering an al- 
ternative aspect of photoelectric heating efficiency. The 
efficiency of the pho toelectric effect on a grain i s tradi- 
tionally defined (e.g. iHollenbach fc Tielensl[T999T ) as the 
ratio of gas heating to the grain UV absorption rate. 
Observationally, this is often approximated as the ra- 
tio of [CII] 158/im or [CII]158/im+[OI]63/xm to the in- 
frared emission from dust, assuming that [CII]158/xm and 
[OI]63/xm account for most of the cooling of the heated 
gas, and total infrared emission accounts for most of the 
UV absorbed by grains. However, photoelectric heating 
is dominated by the smallest grains in the ISM ([Watson! 
Il972t I Jural fl97a ). whereas the total infrared emanates 
prima rily from th e large r grains. To avoid this complica- 
tion, (Helouiran $2QM) used [CII]158/xm/F(5-10/xm) as 
an estimate of efficiency limited to PAH, and showed that 
it varies considerably less than [CII]158/xm/FIR, thus 
demonstrating that the smallest grains, and more specif- 
ically PAHs, do indeed dominate photoelectric heating. 
However, these definitions do not account for the shift- 
ing balance of the diverse dust species in mass ratio or 
absorption of heating radiation. In Figure [13] we see 
a clear difference in (/pah, the mass fraction of PAHs 
to large grains, between the starburst ring and Enuc 
S, much larger than the difference in ([CII]158/iin + 
[OI]63/xm)/PAH. A KS test revealed that the difference 
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in gpAH between the ring and E nuc S is significant. Using 
dust mass and (/ pah maps fromlAniano et all ([2012) , con- 
structed using (|Draine et al, 2007) models, we are able 
to measure the total PAH mass in the ring and Enuc S 
and estimate the amount of [CII]158/im and [OI]63/im 
line emission per unit PAH mass in the ring and Enuc 
S. The values for ([CII]158/im+[OI]63Mm)/M Pj4ff are 
7QLq/Mq for enuc S and 170i©/M Q for the ring, and 
may be interpreted as rough estimates of the amount of 
photoelectric heating contributed per unit mass of PAHs. 
These values are for PDR gas only, i. e., the ionized 
gas fraction of [CIi]158/zm was subtracted before deriv- 
ing these quantities. These ratios suggest that per unit 
mass, PAHs in the ring are responsible for generating 
significantly more, by about a factor of two, [CIi]158/im 
and [OI]63/im emission than PAHs in Enuc S, in spite 
of the slightly lower ([CII]158/im+[OI]63/zm)/PAH flux 
ratio. In other words, the photoelectron production rate 
per unit mass of PAH is higher in the ring, even though 
the photoelectric heating efficiency is slightly lower. This 
difference in ratios is simply a reflection of the greater 
heating intensity available to the PAHs in the ring, mak- 
ing it possible for a PAH grain to process more photons 
per unit time, even if the photoelectric heating efficiency 
is lower, possibly because of increased ionization of the 
PAHs, or because of more frequent two-photon events. 
Since the PAHs are excited stochastically by one photon 
at a time, the increased production rate per unit mass 
can be seen as the result of approximately doubling the 
PAH excitation rate. This reflects how the mean starlight 
intensity (U) va ries between the ring and Enuc S. Using 
(U) maps from lAniano et all (|2012j ). constructed using 
( Draine et aJJ 120071) models, we measured that (U) in 
the ring is ~ 3.6 times higher than in Enuc S, with an 
uncertainty of ~ 30%. Given the uncertainties involved, 
the ratio of (U) and the ratio of photoelectric produc- 
tion rate per unit mass between the ring and Enuc S are 
consistent. 

4.3. What is exciting the H2 ? 

We have seen that the lower [CII]158/im/PAH observed 
in the ring is partially compensated by the increase of the 
[OI]63/im line, as it becomes a more efficient coolant in 
denser and warmer regions. However, another mecha- 
nism that may be respons ible for gas heating in these 
regions is shock heating. iHollenbach fc McKed (1989) 
predicts fast (> 30 kms" 1 ) J shocks to have strong 
[OI] emission. According to these models, the [OI]63/im 
line is typically stronger than the [CII]158/xm line in 
regions where J-shock heating is important. However, 
[OIl63 /nm/[CII1158/xm v aries between 0.3 — 0.5 in the 
ring (|Beirao et al.l 12010). Therefore fast shocks are not 
the dominant heating source of the ISM in the ring, at 
least not at scales probed by the PACS beam, although 
they could still be responsible for much of the [OI]63/im 
emission. Slow (< 30 kms -1 ) C shocks could still be 
resent in the ring, as they emit only weak [OI]63^m 
Hollenbach fc McKed fl989f) . 

By using the radiation field estimates from [CII]158/im 
and [OI]63^m, we can predict the H2 emission from 
PDRs. Assuming the radiation field intensity Go cal- 
culated in section [4~Tl the H2 fluxes from the ring can be 
reproduced using a pure PDR model by Kaufma n et al.l 
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Fig. 13.— Upper: Plot of the ([CII]158/Ltm+[Ol]63/irn)/PAH 
(5.5 — 14fim) ratio as a function of F„(8^tm)/F„ (24^tm) as tr aced by 
Spitzer-YRAC bands for the regions defined in Section l3.3l Lower: 
Plot of the ([CII]158Atm+[OI]63Atm)/PAH (5.5 - 14/nm) ratio as a 
function of the PAH m ass abundanc e qpAH > calculated using the 
dust emission models of Draine k, Li (2013). 



(|2006l e.g.) with high density molecular gas (n ~ 
10 5 cm~ 3 ). Such a mo l ecular gas density is quite plau- 
sible, as iKohno et al.l (|2003| ) reported the existence of 
HCN (1-0) emission in the nucleus and ring, for which 
molecular gas densities of ~ 10 5 cm -3 are necessary (e.g. 
iMeiierink et"al|[200l . For G = 10 2 25 and n H = 10 35 
calculated in Section 4.1, and assuming that the frac- 
tion of TIR flux from PDRs with these conditions is 
responsible for 25% of the total TIR emission, we pre- 
dict that such PDRs will emit a total of H2 flux of 
F(H 2 S(0) - 5(3)) ~ 2.1 x 10~ 16 Win" 2 . Our measured 
total H 2 flux is Fh 2 ~ 1.3 x 10~ 15 Wm~ 2 , about a fac- 
tor of 6 above the modeled value. Therefore, if most of 
the dust emission is coming from low (U) diffuse ISM, 
it is very difficult to account for the observed H2 from 
the ring with intense PDRs. The H 2 lines must predom- 
inately come from some other type of region. 

In the Milky Way there is increasing evidence that 
the diffuse ISM has pockets of warm H2 that radi- 
ate in the H2 rotational lines, where the excited H2 is 
in excess of what would be produced by UV-pumping. 
Both H2 and PAH emission have been observed in ga s 
clouds illuminated by starlight (e.g. lHabart et aLll2004j ). 
Ingal ls et al.l (|2011l ) used Spitzer-IRS to show that high- 
latitude translucent clouds, illuminated by just the lo- 
cal starlight background ((U) ~ 1), may contain small 
amounts of hot, collisionally excited i?2- In these clouds. 
F(H 2 S(2))/F{PAH7.7) ~ 0.002 on average. In the 
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Fig. 14. — Ratio of H2 luminosity summed over the S(0)-S(3) 
lines to 7.7/im PAH luminosity vs. uF v (2A fim) continuum lumi- 
nosity. This ratio indicates the relative importance of mechanical 
heating and star formation power. The open square is the aver- 
age for the nucleus, and the full s quar es are the values for three 
regions in the ring shown in Figure [TOl All o ther points are ratio s 
calculated for galaxies in the SINGS sample (Rousscl ct al. 2007). 

ring of NGC 1097, F(iJ 2 5(2))/F(PAH7.7)~ 0.0015, so 
translucent clouds with low incident radiation are a 
plausible origin for the observed H2 emission in the 
ring. The warm regions that prod uce the H 2 emis- 
sion observed by llngalls et al.l (|2011f ) will also radiate 
in [CH1158/n ni and [01163/xm. Combining the H2 mea- 
sureme nts inllngalls et al. | (1201 If) with the [CII]158/im re- 
sults of llnealls et al.1 (120021) gives F(# 2 5(0)-5(2))/[CII] 
~ 0.4 — 1 for Milky Way translucent cloud positions 
detected in _ff 2 emission, which is much higher than 
the o bserved F(H?.S(0) - 5 (2)) /[CUV 0.090 in the 
ring. In galls fc Dra inc (2012) predict that the translu- 
cent clouds should have [OI] /F(H 2 S(0) - 5(2)) < 0.19, 
well below the observed [OI]/F(H 2 S{0) - 5(2)) - 4.5 in 
the ring. Thus, it seems possible that much of the H 2 
rotational line emission from the ring could come form 
warm regions in the H2 clouds, but these regions do not 
dominate the [CII]158^m or [OI]63/im emission we ob- 
serve. 

In Figure Q3] we plot the ratio of H2 luminosity in 
the 0-0 S(0)-S(3) lines over Lpaht.7 versus F24 for re- 
gions 1, 2 , and 3 in the ring and in the nucleus and 
compare these measuremen ts with data from th e SINGS 
sample of nearby galaxies (|Roussel et al.l l2007h . While 
the values fo r the ring are simila r to those of star form- 
ing galaxies (jRoussel et al.l [2007h . which have a median 
ratio F{H 2 )S{Q) - 5(3)/L(PAH7.7) = 0.011, the nucleus 
is slightly larger, with a value of ~ 0.02, more consis- 
tent w ith the LINER and Seyfert nuclei in Rousscl et al. 
(2007), although there is a large scatter among the galax- 
ies with AGN. 

The nucleus of NGC 1097 has a weak AGN, appar- 
ently fed by an inflow of gas inside the ring ()Prieto et al.l 
2005) . It is possible therefore, that X-rays from the AGN 
might contribute to the heating of the H2 gas in the 
immediate vicinity of the nucleus. The cooling by H2 
rotational lines in X-ray-do minated regions (XD Rs) is 
2% of the total gas cooling ([Malonev et al.| [T996) for a 
temperature of 200 K, typical of warm H2 that emits 
in the mid-infrared. At this temperature, the ratio of 



the first four rotational lines to the total rotational line 
luminosity is L(H 2 S(0) - S(3))/L(H 2 ) = 0.58. Com- 
bining the above factors, we estimate a maximum H2 to 
X-ray luminosity ratio of L(H 2 S(0) - S(3))/L x (2 - 10 
keV)= 0.01. This ratio is conservative, since it assumes 
that all of the X-ray flux from the AGN is absorbed 
by the XDR. The nucleus of NGC 1097 has an X-ray 
luminosity of L x (2 - 10 keV) = 1.73 x 10~ 15 W mT 2 
(|Lutz et al.ll2"00l . Therefore the maximum H2 luminos- 
ity should be L(H 2 ) - 1.7 x 10~ 17 W m~ 2 , significantly 
lower than the measured L(H2S(0)~ 5(3)) in the nucleus, 
which is - 2.3 x 10~ 16 W m~ 2 . Considering L(H 2 S(0)- 
S(3))/L(i7 2 )= 0.58, the contribution of X-rays for heat- 
ing the H2 is less than 5%. This is consistent with the 
near infrared r otational-vibrational sp ectra taken across 
the nucleus bv lReunanen et all (|2002f) . 

Another possible source for H2 heating are cosmic rays 
from supernovae. The energy released by the dissoci- 
ation of each H2 molecule by cosmic rays is ^12 eV, 
including the contrib ution from Ht recom bination and 
H2 re- formation (e.g. iLe Petit et ail 120061 ). To balance 
line cooling, the ionization rate per H2, Ch2, must be few 
x 10~ 14 s _1 . For such a rate, cosmic rays are the main 
destruction path of H 2 molecules. From the Table [5] the 
total H2 emission-line luminosity is 5.91 x 10 40 erg s _1 for 
the ring and 1.02 x 10 40 erg s _1 for the nucleus. Dividing 
this luminosity by the warm H2 masses of 1.52 x 10 s M Q 
and 9.6 x 10 6 M Q respectively, we can estimate a cooling 
rate through the H 2 lines of 3.2 x 10~ 25 erg s -1 H^ 1 
and 8.8 x 10 -25 erg s _1 H^ 1 respectively. The mass 
fraction of the gas in the molecular state depends on 
the rate of ionization and gas density Ch2/"H2- Divid- 
ing the H 2 S(0)-S(3) luminosity by the energy released 
per ionization, we estimate the required ionization rate, 
Ch 2 ~ 1 x 10~ 13 s" 1 for the ring and Ch 2 ~ 6 x 10~ 13 s" 1 
for the nucleus. This value is larger than what is in- 
ferred from H$ observations in the Milky Way, by a 
factor 200 — 1000 fo r the diffuse interstellar medium 
(|Indriolo et al.l 120071 ) . and by more than a factor 20 for 
the molecular gas within 200 pc from the Galactic center 
(|Goto et al.H2008l ). It is therefore unlikely that the warm 
molecular gas in the ring and nucleus of NGC 1097 is 
heated solely by cosmic rays. 

5. SUMMARY 

In this paper we have attempted to understand the 
energetics and the physical conditions in the ISM in 
NGC 1097, using Herschel-PACS and Spitzer-IRS. Our 
goal was to study how heating and cooling of the ISM 
vary with the ISM phases, studying the origin of the 
variations of heating efficiency with dust temperature, 
and the properties of the warm H2 in regions with dif- 
ferent ISM heating mechanisms. We focused on the 
central region, where a ring of star formation and an 
AGN are located, and in two extranuclear regions at 
the ends of the bar, which we named Enuc N and S. 
We divided the ring, Enuc N and Enuc S regions into 
12" x 12" regions. After subtracting the fraction of 
[CII]158/im flux emitted from ionized gas, estimated as 
~ 33% in the ring and ~ 5% in enuc S, we exam- 
ined the correlations between the [CII]158/im/PAH ra- 
tio and the 70/mi/100/im, the PAH ra tio 11.3/7.7/xm. 
We implemented PDR models by Kaufma n et al.l 
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(2006), and used the [CII]158/xm/[OI]63/xm ratio and 
([CII]158/xm+[OI]63/im)/FrR to plot a grid of the lo- 
cal FUV flux Go, and the PDR density n and diagnose 
the properties of the ISM on the ring and Enuc S. We 
also calculated the total H 2 and the PAH 7.7/zm fluxes in 
three regions in the ring and one centered in the nucleus 
to determine the origin of the warm Hi excitation. As a 
result of our analysis we find: 

• Assuming that PDRs account for all [OI]63um, FIR 
flux and 67% of the [CII]158/im, the average inten- 
sity of the radiation field Go is ~ 80 in Enuc S and 
~ 800 in the ring. The gas density is 10 3-3 - 5 cm~ 3 
in both regions. If instead 75% of the FIR flux in 
the ring arises from low starlight intensity regions 
and only 25% of the FIR flux arises from PDRs, 
then Go is ~ 10 23 and tih ~ 10 3 5 in the ring. 

• The ratio [CII]158/mi/PAH(5.5 - Ufxm) is a fac- 
tor of ~ 1.7 lower in the ring than in the Enuc 
regions and [CII]158/im/FIR is a factor of ~ 4 
lower in the ring than at the southern edge of 
the bar. The average 11.3/7.7^m PAH ratio is 
lower in the ring than at the ends of the bar, sug- 
gesting a larger fraction of ionized grains in the 
starburst ring. However, no clear correlation was 
seen between the 11.3/7.7/mi PAH ratio and the 
[CII]158//m/PAH(5.5 — 14/im) ratio across the com- 
plete dataset, implying that grain ionization is not 
the sole driver of the grain heating efficiency vari- 
ations. The ratio ([CII]158^m+[OI]63/im)/PAH is 
still lower in the ring than the ends of the bar, al- 
though the difference is much smaller and the scat- 
ter larger than in [CII]158/mi/PAH(5.5 - 14/zm). 
After the removal of the ionized gas component of 
the [CII]158/im emission and adding [SiII]34.5/im, 
we calculated a gas heating efficiency in the ring 
of ~ 5.2% and ~ 7.7% in Enuc S. The amount of 
[CII]158/im and [OI]63/im line emission per unit 
PAH mass is roughly a factor of two larger in the 
ring than in Enuc S. Since the PAHs are excited 
stochastically by photons, the increased produc- 
tion rate per unit mass can be seen as the result of 
approximately doubling the PAH excitation event 
rate. 

• Using the [CII]158/mi and [OI]63/im emission lines 
and total FIR emission from the ring are used 
in conjunction with standard PDR models, the 
derived radiation field can account for nearly all 
the warm H2 emission measured with Spitzer/IRS. 
However, if we incorporate the results of the dust 



SED modeling, which suggests nearly 75% of the 
FIR flux comes from dust heated by low starlight 
intensities in the diffuse ISM, there is a large excess 
of H2 which cannot be excited in PDRs. In anal- 
ogy to what is seen in the ISM of the Milky Way 
(Inga lls et al.ll20fTI ). we suggest this H2 emission in 
the NGC 1097 ring may be coming from pockets of 
warm H2 gas inside translucent molecular clouds 
in the diffuse ISM. 

• The H 2 S(0)-S(3)/PAH7.7^m ratio is enhanced in 
the nucleus (the site of a weak AGN) relative to the 
ring by a factor of two, as is the fraction of warm 
(T ~ 400 K) molecular gas. X-rays and cosmic 
rays cannot provide more than 1-5% of the energy 
necessary to produce the H2 emission in the nu- 
cleus. Shock waves, possibly associated with the 
central AGN, may be an important contributor to 
gas heating in the nucleus. 

This paper is an example of the study that can be 
carried out by carefully combining Herschel/PACS and 
Spitzer/IRS data focusing on starburst rings or other 
resolved regions of enhanced star-formation in nearby 
galaxies. In future papers we will explore the energy 
balance in the ISM of local star-forming galaxies within 
the KINGFISH sample, and apply these methods of link- 
ing the atomic, and molecular gas feature emission with 
the dust properties, to more luminous starburst galaxies 
observed with Spitzer and Herschel. 
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TABLE 2 

Emission line fluxes from Spitzer-IRS and PACS spectra (in units of 10~ 17 Wm -2 ) 



Region 


Aperture Size 


PAH 


[Aril] 


PAH 


PAH 


[Arlll] 


[SIV] 


PAH 


PAH 


[Nell] 


[Nelll] 


[SIII] [SIII] 


[Sill] 


PAH total 


PAH ratio 


PAH ratio 






6.2 fim 


6.99/^m 


7.7^m 


8.6fim 


8.99/^m 


10.5 M m 


11.3/j.m 


12.6^m 


12.8Mm 


15.6 M m 


18.7Mm 33.5^m 


34. S^m 


6.2 - 12.6/im 


11.3/7.7 


6.2/7.7 




arctscc 2 




15.8 cV 






27.6 cV 


34.8 cV 






21.6 eV 


41.0 cV 


23.3 cV 23.3 cV 


8.2 cV 









laj Nucleus 96.8 2.44 0.103 11.2 1.93 0.275 3.69 2.34 0.364 0-415 0.136 0.246 0.576 21.6 0.33 0.23 

•fcl ±0.19 ±0.025 ±0.6 ±0.12 ■■■ ±0.003 ±0.07 ±0.07 ±0.004 ±0.050 ±0.001 ±0.003 ±0.003 

p2] Ring 680.1 32.6 1.726 130.4 23.3 0.518 0.206 30.5 20.3 4.59 0.205 1.18 1.72 3.81 237.1 0.23 0.25 

±1.3 ±3.68 ±4.6 ±1-0 ±0.013 ±0.008 ±0.5 ±0.7 ±0.01 ±0.003 ±0.01 ±0.01 ±0.01 

Enuc N 1067.2 1.68 0.048 5.79 - - ■ 0.009 1.03 1.32 0.129 ■ ■ • ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ • 10.5 0.23 0.29 

±0.09 ±0.009 ±0.08 ■ ■ ■ ±0.002 ±0.15 ±0.11 ±0.002 
Enuc S 739.3 1.59 0.074 5.26 1.06 ■ ■ ■ ■ ■ ■ 1.32 0.66 0.178 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ 9.89 0.25 0.30 
±0.09 ±0.002 ±0.07 ±0.17 ■_■ ■_■ ±0.10 ±0.06 ±0.35 ■_■ ■_■ ■_■ ■_■ ■_■ 
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TABLE 3 

H2 LINE FLUXES a FROM LOW AND HIGH-RESOLUTION SPECTRA 



Region 


Area (arcsec 2 ) 


S(0)28.2/im 


S(l)17.0/an 


S(2)12.3/mi 


S(3)9.7/im 


Nucleus 


96.8 


0.55±0.13 


6.13±0.78 


6.21±0.79 


10.1±1.01 


Ring 


680.1 


10.7±0.79 


35.9±1.2 


20.7±1.4 


58.5±16.1 


Ring 1 


96.8 


0.77±0.28 


4.00±0.63 


6.13±0.78 


11.7±1.1 


Ring 2 


96.8 


1.96±0.44 


4.26±0.65 


4.53±0.67 


4.66±0.68 


Ring 3 


96.8 


0.99±0.31 


4.44±0.67 


5.86±0.77 


7.17±0.85 


Enuc N b 


1067.2 


0.209±0.039 




0.04±0.01 


0.561±0.010 


Enuc S b 


739.3 


0.322±0.074 






0.770±0.015 



a In units of 10~ 17 Wm" 2 
b From high-resolution spectra 



TABLE 4 

PAH BAND AND LINE FLUXES a FOR EACH 12" REGION 



Regie 




PAH 7.7(im 


PAH 11.3fim 


PAH 6.2 - 12.6^m 


[SIII]18.7Mm b 


[SIII]33.5,.im b 


[SiII]34.8um b 


[OI]63um 


[CII]158um 


FIR C 


Ring 






















1 




177±S 


42.9±0.7 


366±7 


2.02±0.09 


4.15±0.01 


7.19±0.02 


4.38±0.45 


9.91±0.03 


2280 


2 




256±9 


83.9±1.0 


545±5 


2.61±0.10 


4.50±0.01 


8.52±0.02 


7.37±0.21 


17.4±0.1 


3363 


3 




195±6 


50. Oil. 


423±6 


2. 01±0. 11 


3.95±0.01 


7.46±0.02 


6.06±0.18 


14.2±0.1 


2467 


4 




244±7 


59.7±1.1 


514±7 


2.22±0.10 


4.34±0.01 


8.23±0.02 


5.44±0.45 


15.3±0.1 


3174 


5 




235±7 


60.6±1.0 


508±7 


2.30±0.10 


4.00±0.01 


8.30±0.02 


6.64±0.35 


17.4±0.1 


3317 


6 




171±5 


43.7±0.8 


358±6 


1.45±0.14 


3.00±0.01 


6.88±0.01 


3.46±0.27 


12.6±0.1 


2215 


7 




226±8 


57.3±1.1 


474±9 


1.92±0.13 


3.63±0.01 


7.99±0.02 


5.49±0.55 


17.5±0.1 


3101 


8 




172±8 


44.1±0.8 


360±6 


1.68±0.13 


3.05±0.01 


7.14±0.02 


4.42±0.45 


12.8±0.1 


2318 


Enuc 
1 


N 


7.30±0.35 


2.02±0.05 


15.9±0.3 


0.359±0.007 


0.672±0.107 


0.629±0.038 




0.920±0.046 


71.4 


2 




6.73±0.43 


1.92±0.05 


14.7±0.3 


0.291±0.006 


0.642±0.102 


0.601±0.036 




0.712±0.035 


56.9 


3 




7.40±0.38 


2.16±0.06 


15.1±0.3 


0.221±0.005 


0.526±0.084 


0.486±0.015 




0.801±0.040 


75.4 


4 




7.54±0.39 


2.05±0.04 


17.0±0.4 


0.232±0.005 


0.703±0.111 


0.655±0.039 




0.889±0.044 


71.0 


5 




6.71±0.52 


1.73±0.06 


15.7±0.4 


0.117±0.003 


0.573±0.091 


0.535±0.017 




0.749±0.037 


56.9 


6 




9.75±0.28 


2.94±0.05 


22.0±0.3 




0.613±0.018 


0.656±0.104 




1. 01±0. OB 


91.3 


7 




6.99±0.45 


2.05±0.06 


16.2±0.4 




0.541±0.032 


0.582±0.093 




0.819±0.040 


71.2 


8 




8.43±0.56 


3.01±0.06 


18.2±0.3 










0.936±0.042 


103 


9 




8.86±0.56 


2.86±0.06 


21.2±0.4 




0.713±0.113 


0.671±0.040 




1.11±0.05 


97.3 


10 




5.32±0.69 


1.67±0.09 


13.4±0.3 




0.458±0.073 


0.427±0.026 




0.676±0.033 


54.1 


Enuc 
1 


S 


8.96±0.48 


2.65±0.04 


20.7±0.4 








0.360±0.102 


1.12±0.03 


118 


2 




10.4±0.5 


2.82±0.04 


23.7±0.4 




0.068±0.01 


0.315±0.09 


0.356±0.046 


1.31±0.01 


104 


3 




11.6±0.5 


3.18±0.04 


26.5±0.3 








0.512±0.058 


1.86±0.02 


148 


4 




10.3±0.5 


2.60±0.04 


23.0±0.3 


0.058±0.024 


0.071±0.023 


0.229±0.071 


0.232±0.036 


1.10±0.01 


78.2 


5 




13.0±0.6 


3.24±0.04 


28.8±0.3 


0.068±0.032 


0.110±0.032 


0.436±0.131 


0.364±0.036 


1.74±0.01 


122 


6 




5.45±0.51 


1.50±0.05 


12.7±0.5 








0.081±0.013 


0.688±0.007 


53.4 


7 




7.50±0.60 


1.93±0.05 


17.0±0.5 




0.056±0.02 


0.305±0.108 


0.157±0.020 


0.945±0.010 


67.2 


units of 10 -16 Wm" 


2 

















b This line was measured from SH ([SIII]18.7/xm) and LH maps ([SIII]33.5/mi and [SiII]34.8^tm). The SH maps have an area of 1164 
arcsec 2 for the ring, 684 arcsec 2 for Enuc N and 337 arcsec 2 for Enuc S. The LH maps have an area of 1392 arcsec 2 for the ring, 2188 
arcsec 2 for Enuc N and 1193arcsec 2 for Enuc S. 

c The FIR flux was calculated as FIR=TIR/2, with TIR calculated from IDale fc Heloul II2002T ) using Spitzer MIPS 24/^m and Herschel 
PACS 70/im and 160/mi 



TABLE 5 
H2 Model Parameters 



Regions 


Area (arcsec 2 ) 


T(K) 


Ortho/Para 


N{H 2 T 


M(H 2 ) b 


L{H 2 ) C 


Nucleus 


96.8 


100 


1.587 


8.71±3.76 


9.14±3.94 


0.06 






412±19 


3.000 


0.46±0.08 


0.48±0.09 


2.58 


Ring total 


680.1 


119±11 


1.972 


16.0±5.0 


151.0±47.2 


2.75 






467±79 


2.999 


0.15±0.08 


1.44±0.72 


12.6 


Region 1 


96.8 


100 


1.587 


14.1±7.1 


14.8±7.48 


0.09 






470±27 


3.000 


0.29±0.06 


0.30±0.06 


2.73 


Region 2 


96.8 


100 


1.587 


41.3±10.8 


43.4±11.3 


0.28 






375±23 


2.995 


0.40±0.11 


0.42±0.11 


1.52 


Region 3 


96.8 


100 


1.587 


18.0±7.8 


18.9±8.2 


0.12 






401±22 


2.998 


0.41±0.09 


0.43±0.10 


2.05 



a In units of 10 20 cm 2 
b In units of 10 6 M Q 
c In units of 1O 6 L 



